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Along with the vast progress in experimental quantum technologies there is an increasing demand for the 
quantification of entanglement between three or more quantum systems. Theory still does not provide 
adequate tools for this purpose. The objective is, besides the quest for exact results, to develop operational 
methods that allow for efficient entanglement quantification. Here we put forward an analytical approach 
that serves both these goals. We provide a simple procedure to quantify Greenberger-Horne-Zeilinger-type 
multipartite entanglement in arbitrary three-qubit states. For two qubits this method is equivalent to 
Wootters' seminal result for the concurrence. It establishes a close link between entanglement quantification 
and entanglement detection by witnesses, and can be generalised both to higher dimensions and to more 
than three parties. 

I t is a fundamental strength of physics as a science that most of its basic concepts have quantifiability built into 
I their definition. Just think of, e.g., length, time, or electrical current. Their quantifiability allows to measure and 
? compare them in different contexts, and to build mathematical theories with them\ There is no doubt that 
entanglement is a key concept in quantum theory, but it seems to resist in a wondrous way that universal principle 
of quantification. The reason for this is, in the first place, that entanglement comes in many different disguises 
related to its resource character, i.e., what one would like to do with it. In principle, there are numerous task- 
specific entanglement measures^'^. However, most of them cannot be calculated easily (nor measured or esti- 
mated) for generic mixed quantum states, and therefore it is difficult to use them. 

There are notable exceptions, the concurrence^ and the negativity for bipartite systems^. These measures have 
already provided deep insight into the nature of entanglement, but they also have their shortcomings. The 
concurrence is strictly applicable only to two-qubit systems while for the negativities it is not known how to 
distinguish entanglement classes. The generalisations of the concurrence (such as the residual tangle^) do quantify 
task-specific entanglement even for multipartite systems but again it is not known how to estimate them for 
general mixed quantum states. 

There is another difficulty. An AT-qubit density matrix p is characterised mathematically by 2^^ — 1 real 
parameters. Reducing it to its so-called normal form^ — which contains the essential entanglement informa- 
tion — removes 6N parameters. The entanglement measure is determined by the remaining exponentially many 
parameters which need to be processed to calculate the precise value. Even an operational method similar to that 
of Wootters-Uhlmann^'^ would quickly reach its limits with increasing N. Therefore it is desirable to develop 
methods which provide useful approximate answers even for larger systems. If one asks for mere entanglement 
detection, witnesses^^ are such a tool because here the number of required parameters (both for measurement and 
processing) can be reduced substantially. There are also estimates of entanglement measures using witness 
operators^ ^'^^ which, however, have not yet produced practical methods for entanglement quantification. 

Here we develop an easy-to-handle quantitative witness for Greenberger-Horne-Zeilinger (GHZ) entangle- 
ment^^ in arbitrary three-qubit states. It yields the exact three-tangle for the family of GHZ-symmetric states^^, 
and those states which are locally equivalent to them. For all other states, the method gives an optimised lower 
bound to the three-tangle. Due to this feature we call the approach a witness. 

We start by defining the GHZ symmetry^^ and stating our central result. Then we prove the validity of the 
statement for two qubits. We obtain a method equivalent to that of Wootters-Uhlmann, i.e., it gives the exact 
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concurrence for arbitrary density matrices. Subsequently we explain 
the extension of the approach to arbitrary three -qubit states. 

Results 

The procedure. The AT- qubit GHZ state in the computational basis is 

defined as |GHZ> = — (|00 . . . 0) + 1 11 . . . 1)). It is invariant under: 

(i) Qubit permutation, (ii) Simultaneous spin flips i.e., application of 
(jf^. (in) Correlated local z rotations: 

= e^^i (X) e^^^^^ (X) . . . e " ' ' ^^■) ( 1 ) 

where a^, Oy, are Pauli matrices. An AT- qubit state is called GHZ 
symmetric and denoted by if it remains invariant under the 
operations (i)-(iii). An arbitrary AT-qubit state p can be 
symmetrized by the operation 

p\p)=^dUGHzUGHZpUl^^ (2) 

where the integral denotes averaging over the GHZ symmetry group 
including permutations and spin flips. Notably, the GHZ-symmetric 
iV-qubit states form a convex subset of the space of all AT-qubit states. 

Observation: If an appropriate entanglement measure fi is known 
exactly for GHZ-symmetric N-qubit states p^, it can be employed to 
quantify GHZ-type entanglement in arbitrary N-qubit states p. Here, 
^(i/^) is a positive SL{2,C)®^ -invariant function of homogeneous 
degree 2 in the coefficients of a pure quantum state il/, and /nip) is 
its convex-roof extension^^. The estimate for fi(p) is found in the 
following sequence of steps: 

(1) Given a state p, derive a normal form p^^(p), i.e., apply local 
filtering operations so that all local density matrices are proportional 
to the identity'^ {see Section Methods). If p^^{p) = 0 the procedure 
terminates here, and fi(p) = 0. 

(2) Renormalise p^Vtr p^^ and transform it using local unitaries 
V E SU(2)®^ to obtain the state 

NF 

according to appropriate criteria (see below) so that the entanglement 
o/p'(p^7tr p^^) is enhanced. 

(3) Project the state onto the GHZ-symmetric states p^^(p)h^ 
p^(p^^). The estimate for /nip) is obtained after renormalisation 

Two qubits. For two qubits the entanglement measure under 
consideration is the concurrence C(p) (Refs. 4,8). From the 
symmetrization p^{p) of an arbitrary two-qubit state p we flnd (for 
details see Supplementary Information): 

C(p) > max(0, |poo,ii + Pii,oo | + Poo,oo + Pii,ii - 1) • (3) 

In the symmetrization entanglement may be lost, as illustrated by the 

state ) = ( |01) — 1 10) ) for which inequality (3) gives the poor 
v2 

estimate C(^") > 0. Therefore, the optimisation steps (1) and (2) are 
necessary to avoid unwanted entanglement loss in the symmetri- 
zation (3). The goal is to augment the right-hand side of inequality 
(3) up to the point that equality is reached. We will show now that for 
two qubits this can indeed be achieved. 

It is fundamental that the maximum of an SL(2,C)®^-invariant 
function /i(p) under general local operations can be reached by 
applying the optimal transformation pi-^ApA'^/tr ApA"^ where A 
= Ai ® ...®A^ and Aj e SL(2,C) is an invertible local operation^. 
Consider flrst the normal form p^^ip) which is obtained from p 




Figure 1 | The triangle of GHZ-symmetric three-qubit states. The upper 

corners correspond to |GHZ+ ) = ^ (|000) + | 111)) and the lower corner 

to p^(OOl), cf. Ref. 14. The grey area shows GHZ-class states (13 > 0) whereas 
the yeUow area comprises states with vanishing 13 (" W). The border 
between GHZ-class and W-class states is the GHZ/Wline, equation (10) (red 
solid line). We also show a state p^(xoy yo) together with the point (x^, y^) 
that is required to determine the three-tangle T3(:^, 7o)> equation (6). 

by iterating determinant-one local operations^ (see also Methods). 
Such operations (represented by SL(2,C) matrices) describe stochastic 
local operations and classical communication (SLOCC). Consequentiy, 
the normal form is locally equivalent to the original state p, that is, it 
lies in the entanglement class of p. Note that the iteration leading to the 
normal form minimises the trace of the state. Subsequent renormalisa- 
tion increases the absolute values of all matrix elements in equation (3). 
Here, the correct rescaling of the mixed-state entanglement measure is 
crucial. This is why homogeneity degree 2 of fi(\l/) is required^^'^^. 

Hence, transforming p to its normal form increases the moduli of 
Poo,oo» Poo,iiJ Pi 1,00 J Pi 1,11 (and also the concurrence) as much as 
possible for a state that is SLOCC equivalent with p. The sum of 
the off- diagonal matrix elements in equation (3) reaches its max- 
imum if Poo,ii is real and positive. As this can be achieved by a z 
rotation on one qubit we may consider it part of finding the normal 
form and drop the absolute value bars in equation (3). Then, the sum 
of matrix elements equals, up to a factor 1/2, the fidelity of p^^/tr p^^ 

with the Bell state |0+) = — (|00) + 111)). The question is how 
V 2 

large this fidelity may become. 

To find the answer we transform p^^/tr p^^ to a Bell- diagonal form 
using local unitaries (this is always possible^'^^'^^). If then 
Pocfii < |Po?io| apply another SU(2)®^ operation to maximise 
Poou (^^^ Supplementary Information). The result is a Bell-diagonal 
p^^ with maximum real off-diagonal element P^qh (please note that 
p^^ denotes a normalised state, whereas p^^ is not normalised). 
However, Bell-diagonal two-qubit density matrices with this prop- 
erty can be made GHZ symmetric without losing entanglement^ (see 
also Supplementary Information). 

Hence, our optimised symmetrization procedure (l)-(3) leads to 
the exact concurrence for arbitrary two-qubit states p. In passing, we 
have demonstrated that the concurrence is related via C(p) = max(0, 
2f- 1) • trp^^ to the maximum fidelity/ = (0)+ |p^^|a)+ ) thatcanbe 
achieved by applying invertible local operations to p. 

Three qubits. For three qubits, the GHZ-symmetric states are de- 
scribed by two parameters'^ and therefore form a two-dimensional 
submanifold in the space of all three-qubit density matrices. It turns 
out that it has the shape of a flat isosceles triangle, see Fig. 1. A 
convenient parametrisation is 

^(P) = ^(Pooo,iii+Piii,ooo) (4) 
y{p) = (^Pooo,ooo + Piii,iii - ^ (5) 
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Figure 2 | Illustration of the procedure for finding the three-tangle of a 
general mixed three-qubit state p. In the xy plane, there is the triangle of 
GHZ-symmetric states while on the vertical axis, the three-tangle for each 
GHZ-symmetric state (cf. equation (6)) is shown. Simple projection 
► generates a non-optimal GHZ-symmetric state. The optimisation 



steps (2) move the symmetrization image to p^P^'^ e 
enhanced three-tangle. 



as it makes the Hilbert-Schmidt metric in the space of density ma- 
trices conincide with the Euclidean metric. This way geometrical 
intuition can be applied to understand the properties of this set of 
states. All entanglement- related properties of GHZ-symmetric states 
are symmetric under sign change x <r^ —x as this is achieved by 
applying cr^ to one of the qubits. 

The GHZ-class entanglement of three-qubit states is quantified 
by the three- tangle T3 (Refs. 6,17, see also Methods). For GHZ- 
symmetric three-qubit states p^Xq, yo) the exact solution for the 
three-tangle^° (see also Methods) is 



for Xq <x^ and yo <y^ 



V3 



otherwise 



(6) 



w 

-yo 



where Xq^O and (x^, y^) are the coordinates of the intersection of 
the GHZ/W line with the direction that contains both GHZ+ and 
P^^o^yo) (cf. Fig. 1). The grey surfaces in Fig. 2 illustrate this solution. 

Now we turn to constructing a quantitative witness for the three - 
tangle of arbitrary three-qubit states by using the solution in equation 
(6). As before, the main idea is that an arbitrary state can be symme- 
trized according to equation (2) and thus is projected into the GHZ- 
symmetric states. Again, we assume Pooo,iii real and nonnegative, so 
that x(p) > 0. From Figs. 1 and 2 it appears evident that the entan- 
glement of the symmetrization image p^{p) can be improved by 
moving its point (x(p), y{p)) closer to GHZ+. More precisely, the 
entanglement measure is enhanced upon increasing one of the coor- 
dinates without decreasing the other (cf. equations (3) and (6)). 

In this spirit, finding the normal form in step (1 ) is appropriate as it 
yields the largest possible three-tangle for a state p^^/tr p^^ locally 
equivalent to the original p (cf. Ref. 7). As the normal form is unique 
only up to local unitaries it does not automatically give the state with 
minimum entanglement loss in the symmetrization. Therefore, the 
unitary optimisation step (2) is required to generate the best coordi- 
nates. 

In the symmetrization the information contained in various 
matrix elements is lost. For two qubits, however, the concurrence 
of the optimised Bell-diagonal states depends only on P^q^qq + P^qh 
and the loss of p^iiQ in the symmetrization does not harm. In con- 
trast, the three-qubit normal form depends on 45 parameters. We 



may not expect that T^ip) depends only on two of them and, hence, 
entanglement loss in the symmetrization (3) is inevitable (cf. 
Supplementary Information). Consequently, steps (l)-(3) lead to a 
lower bound for the three-tangle that coincides with the exact T^ip) at 
least for those states which are locally equivalent to a GHZ-symmet- 
ric state. The most straightforward optimisation criterion in step (2) 
is to maximise ^(p^(p^^)). Alternative criteria which generally do 
not give the best T^ip) but can be handled more easily (possibly 
analytically) are maximum fidelity (GHZ+ |p^(p^^)|GHZ+), min- 
imum Hilbert-Schmidt distance of p^(p^^) from GHZ+, or max- 
imum Re p^^Q J J . 

Discussion 

Evidently this approach can be generalised. Therefore we conclude 
with a discussion of some of its universal features. The essential 
ingredients are an exact solution of the entanglement measure for 
a sufficiently general family of states with suitable symmetry, and the 
entanglement optimisation for a given arbitrary state p via general 
local operations. The former determines the border where the entan- 
glement vanishes. The latter ensures an appropriate fidelity of the 
image p^(p) with the maximally entangled state. This reveals a 
remarkable relation between entanglement quantification through 
SL(2, C) invariants and the standard entanglement witnesses which 
we briefly explain in the following. 

A well-known witness for two-qubit entanglement is 



1 



I4 — + ) |. It detects the entanglement of an arbitrary 



W2 



normalised two-qubit state p^"^^ if 



0>tr(p2^^W2 



-(a)+|p2^^|a)+). 



On the other hand, from our concurrence result 

C(p2qb) ^maxTo,^ max^ [2(0)+ |Ap2q^A"^|0+) -tr(Ap2qU'^)] 



>2(a)+|p2q^|0+)-tr(p2q^) 
= -2tr(p2^^W2) 



(7) 



we see, by dropping the optimisation over SLOCC A = Aj ® A2, that 
W2 = — 2W2 is a (non-optimised) quantitative witness for two- 
qubit entanglement. In other words, Wi yields one of the many 
possible lower bounds to the exact result. Analogously it is straight- 
forward to establish the relation between the standard GHZ witness 
3 

W3 = -Is — |GHZ+)(GHZ+ I and the non-optimal quantitative 

witness W3 = — 4W3. The latter represents a linear lower bound 
to the three-tangle obtained via the optimisation steps (l)-(3) (see 
Supplementary Information). 

Finally we mention that our approach can be used without opti- 
misation, i.e., either without step (1 ), or (2), or both. This renders the 
witness less reliable but more efficient. At best it requires only four 
matrix elements (for any N) . We note that, if we apply the witness to a 
tomography outcome the measurement effort can be reduced by 
using the prior knowledge of the state and choosing the local mea- 
surement directions such that the fidelity with the expected GHZ 
state is measured directly. This implements optimisation step (2) 
right in the measurement. 



Methods 

Normal form of an N-qubit state. The normal form of a multipartite quantum state 
is a fundamental concept that was introduced by Verstraete et aU. It applies to 
arbitrary (finite- dimensional) multi-qudit states. Here we focus on AT-qubit states 
only. 

In the normal form of an AT-qubit state p, all local density matrices are proportional 
to the identity. Therefore the normal form is unique up to local unitaries. Remarkably, 
the normal form can be obtained by applying an appropriate local filtering operation 
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p^^ = (Ai (X) . . . ®An)p{Ai (X) . . . (x)Aw)^ 

where Aj e SL(2,C). Therefore p^^ is locally equivalent to the original state p. The 
normal form p^^ is peculiar since it has the minimal norm of all states in the orbit of p 
generated by local filtering operations. Practically, the normal form can be found by a 
simple iteration procedure described in Ref 7. It is worth noticing that GHZ- 
symmetric states - which play a central role in our discussion - are naturally given in 
their normal form. 

Three-tangle of three-qubit GHZ-symmetric states. The pure- state entanglement 
monotone that needs to be considered for three-qubit states is the three-tangle 13(1/^), 
i.e., the square root of the residual tangle introduced by Coffman et aU: 



(8) 



T3((A) =2^1^^1-2^2+4^31, 

='Aooo'Aui +«AooiV^iio + 'Aoio^ioi +^oii^?oo 
^^2 = 'Aooo'Aooi'Aiio^iii +'Aooo^oio^ioi'Aiii + 

+ 'Aooo'Aoii'Aioo'Aiii+'Aooi'Aoio'Aioi'Aiio + 

+ «AoOl'Aoil'AlOO'AllO+^01o'Aoil'AlOO'Al01 

^^3 = 'Aooo'Aiio^ioi^oii + 'Aloo'Aoio'Aool ^n\- 

Here ij/jj^i with j,A;,/ e {0,1} are the components of a pure three-qubit state in the 
computational basis. The three-tangle becomes an entanglement measure also for 

mixed states p= Ylijpi\^j)(^j \ the convex- roof extension^^ 



(9) 



i.e.y the minimum average three-tangle taken over all possible pure-state 
decompositions {pp xjjj]. In general it is difficult to carry out the minimisation 
procedure in equation (9), but there exist various approaches for special families of 
states^^'^^'^°"^^ For GHZ-symmetric three-qubit states, the convex roof of the three- 
tangle can be calculated exactly (see equation (6)). This solution is shown in Fig. 2 and 
can be understood as follows. The border between the W and the GHZ states is the 
GHZ/ W line which has the parametrised form^'* 

8(4- v2)' ^ 4 4-v2 ^'"^ 
with — 1 < V < 1. The solution for the convex roof is obtained by connecting each 
point of the GHZ/ W line (x^,;/^, T3 = 0) with the closest of the points (xghz+ ^ i 2' 

7GHZ+ = — > T3 = 1). That is, the three-tangle is nothing but a linear interpolation 

between the points of the border between GHZ and W states, and the maximally 
entangled states GHZ+. 
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